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Using a single particle density distribution for a system of self-gravitating particles which ul- 
timately forms a black hole, we from a condensed matter point of view derive the Schwarzschild 
radius and by including the quantum mechanical exchange energy we find a small correction to the 
Schwarzschild radius, which we designate as the skin of the black hole. 
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I. INTRODUCTION 



In the evolution of a star, when the radius of a star be- 
comes less than a certain limit, called the Schwarzschild 
radius Rs c h — 2GM/c 2 , where M is the mass of the ob- 
ject; G, Newton's universal gravitational constant, and c, 
velocity of light; neither light nor material particles can 
escape from the star. Thus, one finds from above, the 
Schwarzschild radius corresponds to the situation when 
the escape velocity is equal to the velocity of light, a par- 
ticle (or a photon) coming from a large distance when 
passes near by the black holei^, it is not only attracted 
towards it but also its orbit diverges from a straight 
line. It may so happen that if the particle goes too close 
to a black hole it is likely to be trapped and hence, it can- 
not escape to infinity. As a special case, if the particle 
travels straight to the center of the black hole, it falls in- 
side it and is lost forever. An interesting problem that is 
associated with the formation of a black hole is the final 
collapse of a massive star. This happens when the nuclear 
fuel inside the central core of the star gets exhausted. At 
this stage, it is the dominance of the gravitational attrac- 
tion among the particles within a star over the internal 
outward pressure which makes the star to collapse. A 
black hole, though may be formed from baryons, leptons 
etc; the exterior observer cannot have access to the de- 
tails of the inside of the black hole. The observer probes 
the black hole mass M, electromagnetic charge Q and an- 
gular momentum J. This is refered to as baldness of the 
black hole or as Wheeler describes, a black hole has no 
hair—. But actually it has only three hairs, M, Q, J. We 
in this paper derive the Schwarzschild radius of a system 
of gravitating particles quantum mechanically and show 
that when quantum exchange correction is taken into ac- 
count, there is a thin correction to this radius, which 
we designate as the skin of the black hole. It has been 
shown previously^ that using the quantum field theoretic 
method, by including a single-closed-loop in the self en- 
ergy, a quantum correction to the classical Schwarzschild 
solution of the order of ~ G 2 can be found. This comes 
from the gravity sector. The correction that we find, 
comes from the exchange part of the matter-energy sec- 
tor of the black hole. Our correction is also of the order 
of ~G 2 . 



We here consider a black hole with mass only, a 
Schwarzchildl type black hole. As described above, since 
the black hole forms out of particles due to immense grav- 
itational attraction, we consider the system of particles 
as a many particle system that ultimately forms a black 
hole. We have succeeded in developing a quantum me- 
chanical approach^ for a system of self-gravitating par- 
ticles using Newtonian gravity. We calculate the energy 
content of the system including exchange interaction. In 
order to calculate the total energy of a star, we choose a 
trial single-particle density to account for the distribution 
of particles within star. The form of our single-particle 
density is such that it has a singularity at the origin. Ap- 
plying it to the case of a neutron star, we not only arrive 
at a compact expression for the radius of the neutron star, 
but also obtain an expression for the binding energy of 
the star which varies with the particle number— as N 7 / 3 , 
where N is the particle number. Such a dependence with 
N is in agreement with those of the earlier workers^ The 
aim of the present work is to give a derivation of the so- 
called Schwarzschild radius even without using GTR and 
relativistic quantum mechanics. By accounting for the 
exchange effects due to the interparticle correlations to 
the total ground state energy of the system, we find a 
quantum correction to the Schwarzschild radius. 



II. MATHEMATICAL FORMULATION 

In order to describe a system of N self-gravitating par- 
ticles in absence of any source for radiation, we use a 
Hamiltonian of the form: 

1=1 1=1 j=l,ijij 

where v(\Xi — Xj\) — —g 2 j\X; L — Xj\, is the interparti- 
cle interaction between pair of gravitating particles and 
g 2 = Gm 2 , m being the mass of particle and G be- 
ing the universal gravitational constant. In the present 
case we confine ourselves to the system of neutrons only. 
Since the wavefunction of a neutron star is not known, 
we proceed^ to evaluate the total kinetic energy of the 
system using a model density distribution function for 
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particles (neutrons) within the neutron star. The parti- 
cles being fermions, we use the Thomas- Fermi formula for 
calculating the total kinetic energy of the system. For an 
infinite many-fermion system, the average particle den- 
sity and the fermi momentum of a particle are related to 
each other as: 



n = 



3^ 



(2) 



For a finite system like the star, since the density distri- 
bution is a function of radius vector r, the fermi energy 
of a particle is supposed to be dependent on r. For an 
infinite many-fermion system, the total kinetic energy of 
the system^ is given as 



3 3 h 2 

< KE > inf = -ne F = -rt(37r 2 n) 2/3 — 
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(3^ 2 ) 2 / 3 n 5 / 3 , 



(3) 



In analogy with Eq.Q, for a finite system^, we write, 

3h 2 



< KE >-- 



10m 



(3tt 2 ) 2 / 3 / dr[p(f)} 5 / 3 , 



(4) 



The total energy E of the system is given as 

E =< H > = < KE > + < PE >, (5) 

where 



< pe >= - g 4 1 

2 ./ f-r" 



(6) 



The expression given in Eq.((S]) is written in the Hartree 
approximation. In order to find E, we choose a trial 
single-particle density of the form 



p(f) 



esp[-(5) 1/2 ] 

( r )3 /2 



(7) 



where A is the normalization constant, which is deter- 
mined using the relation 



p(r)dr = N. 



(8) 



As one can see from Eq.©, p(r) is singular at r — 0. In 
general one could choose a single particle density of the 
form 



p(r) = A 



exp[-(jY 



(9) 



where v = 1,2, 3, 4. ..or \>---- Integer values of v 
are not permissible because they make the normalization 



found to be most appropriate, because, as we shall see 
later, it gives the expected upper limit for the critical 
mass of a neutron star—, beyond which black hole forma- 
tion takes place. Any other value of v would give rise 



to a different value for the critical mass. Also because 
if v goes to zero (like 1/n, n — > oo), p(r) would tend 
to the case of a constant density as found in an infinite 
many-fermion system. Having accepted the value v = i , 
the parameter A associated with p(r) is determined after 
minimizing E(X) =< H > with respect to A. This is 
how, we are able to find the total energy of the system 
corresponding to its lowest energy state. After evaluating 
the integral shown in Eq.(|4]) and Eq.©, we obtain 

m 257T 16 \ z 16 A 

Differentiating this with respect to A and then equating it 
with zero, we obtain the value of A at which the minimum 
occurs. This is found as: 



25m 5 2V 16 ; N 1 / 3 



(11) 



Here we are only concerned with the total kinetic energy 
of the system. At A = Ao, we have, 

4 

< KE >= 0.015441^-iV 7/3 (12) 
nr 

The total energy E of the system is found to be just 
negative of this. 



III. DERIVATION OF SCHWARZSCHILD 
RADIUS 

We now try to calculate the average velocity of a parti- 
cle within the neutron star using Eq.(fl"2"|). Let us denote 
it by < v 2 >. If M denotes the total mass of the neutron 
star, one writes 



< KE >= ~M < v 2 >, 



(13) 



where M — Nra. By comparing Eq. (fT3"l) with Eq. (fT2")) . 

one obtains, 



a 

< v 2 >= 0.030882^- N 4/3 . 



h 2 



(14) 



From the expression for the total kinetic energy of an in- 
finite many-fermion system^, one finds that the average 
velocity of a particle within the system is ~ 0.77u/, Vf be- 
ing the fermi velocity^ of the particle within the system, 
which is maximum velocity of that particle. From this, 
one clearly sees that the maximum velocity of a particle 
belonging to an infinite many-fermion system is greater 
than the average particle velocity. In view of this fact, 
we could write the maximum velocity of a particle within 
a neutron star as 

< v 2 > ma x= a < v 2 >= 0.030882a|^A 4 / 3 (15) 

nr 
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where a is a constant whose value is to be greater than 
unity and it is to be calculated later. v max can be identi- 
fied as the escape velocity of a particle within a neutron 
star. According to special theory of relativity, v 2 ^^ is to 
be less than c 2 , c being the velocity of light. That is, 



0.030882a|^iV 4/3 < c 2 , 

h 



From this it follows that 



, T 13. 574409, ftc N o /2 
N ^ 374— ("2 7 



N c (say), 



(16) 



(17) 



having g 2 = Gm 2 . Substituting Eq.(fl7| in Eq.([TT|). one 
finds that, 

Ao > Ac = ^a 1 / 4 [0.8483718(^) 3 / 2 ]. (18) 
Cr g z 

If we define the radius of a neutron star as Rq — 2Ao, we 
have the expression for the critical radius as, 



Rc = 2A C 



2 ^ a l/4 [0 . 8483718( ^ ) 3/2 j 



2GM 



(19) 



Our identification about the radius R of the star with 
2Ao is based on the use of socalled quantum mechani- 
cal tunneling^ effect. Classically, it is well known that a 
particle has a turning point where the potential energy 
becomes equal to the total energy. Since the kinetic en- 
ergy and therefore the velocity are equal to zero at such 
a point, the classical particle is expected to be turned 
around or reflected by the potential barrier. From the 
present theory it is seen that the turning point occurs at 
a distance R = 2Ao- This is the reason why we identify 
2Ao with the radius of a star. For R > 2Ao, a particle, 
belonging to the system, may have an access to the re- 
gion beyond R > 2Ao, because of quantum mechanical 
tunneling, but is forbidden by classical theory. 

R c as given in Eq. (TlT)]) is being identified as the so 
called Schwarzschild radius which we have derived here 
by treating the system as a quantum many-body system. 
When Rq < R Cl the corresponding neutron star becomes 
a black hole. From Eq. fTTj) . we therefore find that the 
lowest mass of the neutron star beyond which black hole 
formation takes place is given as 



M r = mNr 



13.574409 ,hc. 3/2 



,3/4 



r 



(20) 



In order to determine a, we now try to evaluate the lim- 
iting mass of a neutron star following the general expres- 
sion for the radius of a star. Beyond this mass, the black 
hole formation is likely to take place. For that, we con- 
sider the situation when 



(Rq — 2A ) — (R S ch — — 5— ) 
where M — Nm. From this, we arrive at 



(21) 



Since the expression in the right hand side of Eq. ([22|) 
should be equal to the one given in right hand side of 
Eq. (fT7|) . we must have a = 16. Under this situation, we 
have 



■j 2 max = 0.494112^iV c 4 / 3 , 



(23) 



where N c = 1. 696758 (|$) 3 ^ 2 , which, when evaluated, be- 
comes 3.7390777 x 10 57 . For a neutron star in which the 
number of neutrons exceeds N c , it has the tendency of 
forming a black hole. In that case, its mass must exceed 
M = M c = mN c = 3.12213 M©, M Q being the solar 
mass. 



IV. QUANTUM CORRECTION 

So far we have been discussing about the quantum 
mechanical derivation of the Schwarzschild radius Rsch- 
The very form of Rs c h shows that it is a classical re- 
sult, leaving aside the fact the number of particles N 
within a neutron star— is to be less than iV c where 
N c = 1.70(^5-) 3 / 2 , which involves the Planck's constant 
h. Now, inorder to account for the quantum corrections 
to Rsch, we go beyond the Hartree contribution to the to- 
tal energy of the system. That is the exchange correction 
or Hartree-Fock(HF) termii over the Hartree result (di- 
rect contribution). Since the HF-correction term is non- 
local we make use of the local density approximation 11 
to write it as, 



< PE > ex = — (3tO 

Z7T 



2^1/3 g 2 



dr[p{f)} 4/3 . (24) 



This when evaluated gives 



< PE > ex = !(J-)4/3 ( 3, 2) i/3 ff2 ^! (25) 

With the inclusion of this extra term, the expression for 
-E(A), Eq. (|10p is minimized with respect to A and we 
arrive at 

Following the argument discussed earlier, we identify the 
radius of the neutron star by R' = 2X' . As before writing 
v 'max = 16 < v' > and using the condition that v' max < 
c , we obtain 

N < N' c — 1.696758(^) 3/2 [1 + 0.4747761(^)] (27) 



he 



Corresponding to N' c , the new expression for the critical 
radius i?' becomes 



N > (1.696758)( 



he 



3/2 



(22) 



R' = 2A' = 2^^ [1 + 0.7912723(|-)] (28) 
c l he 
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where M c = mN c = 1.696758(|f ) 3/2 . The above expres- 
sion, Eq. (f!?5|) is obtained by keeping terms upto order 
only in Eq. ||2S|). For N > N' c , the neutron star 
is likely to go over the black hole stage. From Eq. ([28|) . 
we find that the second term within the square bracket, 
forms the quantum correction to the Schwarzschild ra- 
dius. As expected, it involves the gravitational fine struc- 
ture constant (j-)- Since it is of the order 10~ 39 , obvi- 
ously it makes an extremely small correction to Rs ch- 
it has been shown earlier— that using the quantum field 
theoretic method and by including a single-closed-loop 
in the self energy, a quantum correction to the classical 
Schwarzschild solution of the order of ~ G 2 can be found. 
This comes from the gravity sector. The correction that 
we get is also of the order ~ G 2 but it comes from the 
exchange part of the matter-energy sector of the black 



hole. 



V. CONCLUSION 

We in this paper derive the Schwarzschild radius of 
a black hole from a condensed matter point of view by 
using a single particle density distribution for the many- 
body self-gravitating system which ultimately forms a 
black hole. By incorporating the quantum mechanical 
exchange interaction, we also find a thin correction to 
the Schwarzschild radius which we designate as the skin 
of the black hole. 

We thank F. Hehl for critically reading the manuscript 
and bringing to our notice the Ref.3. 
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